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1.  SUMMARY 


Broadband  seismic  data  acquired  during  the  Hi-CLIMB  experiment  are  used  to  study  seismic 
events  and  path  propagation  in  the  Nepal  Himalaya  and  the  south-central  Tibetan  Plateau.  The 
2002-2005  experiment  consisted  of  233  stations  along  a  dense  800  km  linear  north-south  array 
extending  from  the  Himalayan  foreland  into  the  central  Tibetan  Plateau.  The  main  array  was 
flanked  by  a  350  km  x  350  km  sub-array  in  southern  Tibet  and  central  and  eastern  Nepal.  The 
dataset  provides  an  opportunity  to  obtain  seismic  event  locations  for  ground  truth  evaluation,  and 
to  detennine  source  parameters. 

The  Hi-CLIMB  experiment  revealed  intense  and  distributed  seismicity  in  the  south-central 
Tibetan  Plateau.  The  Hi-CLIMB  seismic  event  database  contains  over  22,500  events  from  June 
2004  to  August  2005.  Automated  arrival  time  picks  from  vertical  component  data  consist 
primarily  of  P  waves  and  relatively  few  S  waves.  To  increase  the  number  of  S  picks  and  to 
obtain  accurate  picks  at  distances  larger  than  l°-2°,  we  implemented  an  automatic  S-wave  picker 
that  utilizes  signal  attributes  from  three-component  seismic  data.  We  adjusted  picking  parameters 
and  thresholds  to  optimize  automatic  picks  relative  to  manual  picks  from  events  distributed 
throughout  south-central  Tibet.  For  shallow  events,  Sg  can  be  picked  accurately  to  the  Sg/Sn 
crossover  distance  of  about  3°  while  Sn  arrivals  are  absent.  Deep  events  beneath  the  southern 
Tibetan  Plateau  and  the  High  Himalayas  produce  clear  S  arrivals  that  can  be  picked  to  5°-6° 
distance.  We  applied  the  S-picker  to  5790  ML(P)>2.0  events  resulting  in  about  93,000  high 
quality  S  picks. 

We  identified  54  GT5-candidate  events  with  ML(P)>=2.0  of  which  50  satisfy  the  GT5 
criteria  at  the  95%  confidence  level.  Waveform  modeling  resulted  in  143  earthquake  moment 
tensor  solutions. 

2.  INTRODUCTION 

The  broadband  dataset  from  the  2002-2005  Hi-CLIMB  seismic  array  in  the  Himalayan- 
Tibetan  collision  zone  is  unique  in  its  large  aperture  and  dense  station  spacing.  The  purpose  of 
our  research  was  to  improve  source  and  path  calibration  in  regions  of  poor  crustal  propagation  to 
enhance  capabilities  to  estimate  magnitude  and  yield  of  seismic  events  in  low  Q  and  highly 
scattering  environments.  We  have  three  linked  objectives.  The  first  is  to  obtain  high-precision 
ground-truth  seismic  event  hypocenters.  The  second  is  to  provide  seismic  moment  and 
independent  depth  estimates  from  waveform  modeling  of  local  and  regional  earthquakes.  The 
third  is  to  determine  the  distance  evolution  of  seismic  coda,  utilizing  the  quasi-continuous 
control  offered  by  the  dense,  laterally  large  network  (paying  emphasis  on  mantle  coda), 
important  for  low  Q  areas. 

3.  THE  HI-CLIMB  DATASET 

The  Hi-CLIMB  broadband  seismic  experiment  (Nabelek,  et  ah,  2005;  2009)  consisted  of  233 
sites  in  Nepal  and  the  south-central  Tibetan  Plateau  (Figure  1).  The  network  operated  from  2002 
through  2005  and  each  site  was  occupied  for  12  to  20  months.  The  network  consisted  of  an  800- 
km  long  north-south  array  covering  the  India-Eurasia  collision  zone  from  the  Ganges  foreland 
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across  the  High  Himalayas  into  the  Tibetan  Plateau.  Station  spacing  was  3-4  km  in  the  south  and 
about  8  km  north  of  the  Yarlung-Tsangpo  suture.  Complementary  lateral  deployments,  with  30 
to  40  km  station  spacing  in  central  and  eastern  Nepal  and  southern  Tibet  significantly  improve 
seismic  event  location  capabilities,  provide  wave  propagation  control  for  the  transition  from  the 
foreland  into  the  plateau  and  allow  high-resolution  3-D  wave  speed  and  attenuation  tomography. 
The  Hi-CLIMB  data,  recorded  continuously  at  40  and  50  sps,  provide  a  unique  opportunity  to 
detennine  locations,  source  parameters  and  coda  distance  evolution  of  small  sized  seismic  events 
at  unprecedented  accuracy  in  central  Asia. 

Figure  1  shows  the  large  spatial  extent  and  small  station  spacing  of  the  Hi-CLIMB  array. 

This  study  focuses  on  data  from  Phase  2  of  the  project  (red  symbols),  in  an  area  which  is 
relatively  less  known  compared  to  the  Himalayan  region  in  Nepal,  for  which  regional  network 
data  are  available. 


80"  81*  82*  83*  84*  85*  86*  87*  88*  89*  90*  91*  92* 


Figure  1.  Map  of  the  Hi-CLIMB  seismic  network  (phase-2:  red  triangles,  phase-1:  open 
triangles)  and  permanent  GSN  station  Lhasa.  During  the  first  20  months  of  Hi-CLIMB 
(Nepal  and  South  Tibet),  75  broadband  instruments  were  available  while  during  the  second  15 
months  (South  and  Central  Tibet),  111  instruments  were  available.  Altogether,  233  sites  were 
occupied. 
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4.  SEISMIC  EVENT  LOCATION 


4.1.  Event  Detection  and  Initial  Location 

The  Hi-CLIMB  waveform  database  (1.4  TB)  has  been  organized  using  the  Antelope  software 
developed  by  Boulder  Real  Time  Technologies  (BRTT).  The  package  includes  modules  for 
automatic  detection  and  arrival  time  picking,  event-arrival  association,  and  event  location.  Initial 
analysis  of  event  detections  revealed  about  a  thousand  legitimate  event  detections  per  week 
(Figure  2).  Using  a  low  detection  threshold  (STA/LTA  ratio  of  3.0  in  a  frequency  band  of  2-10 
Hz)  and  an  association  rule  of  7  or  more  phases  needed  to  constitute  an  event  resulted  in  an 
initial  database  of  almost  1 10,000  events  for  the  15-month  June  2004  to  August  2005  period. 

This  corresponds  to  an  average  of  240  events  per  day  and  is  more  than  500  times  the  number  of 
seismic  events  listed  in  the  United  States  Geological  Survey  (USGS)  Preliminary  Determination 
of  Epicenters  (PDE)  catalog  for  the  same  region  and  period.  Many  of  these  events  are 
aftershocks  of  the  2004/193  M=6.2  earthquake  on  the  west  side  of  the  array. 


Figure  2.  Event  detection  for  a  20-minute  Hi-CLIMB  data  window  (Julian  day  193,  2004, 
start  at  23:33  UTC).  P-  and  S-phase  association  has  been  performed  and  four  stronger  local 
events  were  recorded  across  the  array  and  many  smaller  events  at  specific  station  groups  only. 
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The  large  number  of  events  rendered  manual  inspection  unfeasible.  We  developed  post¬ 
processing  filters  to  remove  spurious  detections  and  erroneous  phase  and  event  associations  to 
automate  data  management.  Figure  3  shows  a  simplified  flowchart  of  post-processing  steps  that 
implement  simple  logic  considerations  about  what  constitutes  a  valid  seismic  event.  The  main 
considerations  are  (1)  stations  close  to  each  other  have  to  record  an  event  and  (2)  stations  close 
to  the  epicenter  have  to  record  an  event.  The  initial  arrival  set  for  a  declared  event  could  contain 
a  significant  number  of  wrongly  associated  detections  resulting  in  a  grossly  incorrect  location. 
Problems,  for  example,  stemmed  from  wrong  phase  association  caused  by  simultaneously 
occurring  events  in  different  locations  within  the  Hi-CLIMB  array  region,  misplaced  seismicity 
from  the  great  2004  Sumatra  earthquake  sequence  or  incorrectly  associating  P  arrivals  as  S 
arrivals,  or  vice  versa.  The  initial  filter  eliminates  spurious  detections  or  entire  events  based  on 
(1).  In  later  stages,  we  expect  locations  to  be  basically  correct  and  thus  require  (2). 


Figure  3.  Simplified  flowchart  of  database  post-processing  steps  performed  to  clean  up  the 
initial  associated  event  table.  Processing  steps  fall  into  two  groups:  (1)  removing  or 
reassigning  phase  picks  to  obtain  a  more  consistent  set  of phase  picks  for  an  event  leading  to 
improved  location  after  a  relocation  (green  box)  or  (2)  removing  spurious  events  for  which  the 
set  of phase  picks  could  not  be  fixed  (blue  box  “Remove  events”).  Details  are  described  in 
Carpenter,  2010. 

We  manually  picked  P-  and  S-arrivals  to  refine  the  ID  velocity  model  for  location  and  to 
evaluate  automatic  pick  quality.  We  used  two  events  from  the  Pumqu-Xianza  Rift  (PXR)  with 
clear  P  and  S  onsets  across  the  Hi-CLIMB  array  to  refine  the  model  because  their  hypocenters 
are  well  constrained  and  event-station  paths  sample  southern  Tibet  resulting  in  a  representative 
ID  model.  Figure  4  shows  the  observed  and  predicted  (tib3p)  travel  times  for  the  refined  (tib3p) 
ID  model.  The  observed  P  arrivals  are  consistent  with  a  simple  two-crustal-layer  over  mantle 
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model;  the  shallow  low-velocity  layer  agrees  with  receiver  function  analysis  (Nabelek,  et  ah, 
2009)  and  waveform  modeling  (Baur,  2007).  A  single  crustal  S-velocity  layer  is  required  by  the 
traveltime  data  (Figure  4)  since  Sn  crossover  distance  is  beyond  the  Pn  crossover  distance  (at 
about  300  km).  The  upper  mantle  velocities  in  tib3p  are  based  on  the  arrival  times  from  the 
additional  near-Moho  events  from  the  south.  The  model  tib3p,  which  we  unifonnly  used  for 
locating  the  events,  is  strictly  applicable  only  for  the  Lhasa  block.  South  of  the  Yarlung-Tsangpo 
suture  and  north  of  the  Bangong  suture  the  velocity  structure  is  clearly  different,  particularly  the 
depth  of  Moho. 


tibct3  and  tib3p  velocity  models 
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Figure  4.  Right.  Observed  (dots)  and  predicted  (green  lines)  travel  times  vs.  distance  from 
the  epicenter  for  a  manually  picked  PXR  event. 

Figure  4.  Left.  Refined  velocity  model  tib3p  based  on  these  observations  was  used  for  final 
location  of  the  Hi-CLIMB  dataset.  Earlier  analyses  used  model  tibet3,  which  was  estimated 
from  several  type  events  within  study  region.  Red  open  circles  are  analyst-picked  S-arrivals. 

To  assess  the  quality  of  the  automatic  arrival-time  picks  and  to  evaluate  whether  a  single  ID 
model  suffices  for  our  location  purposes,  we  manually  picked  P-  and  S-arrivals  for  16  stronger 
events  distributed  throughout  the  study  region.  Figure  5  shows  differences  between  manual  and 
automatic  picks  for  P  and  S  arrivals.  At  distances  less  than  300  km  for  P  wave  and  less  than  100 
km  for  S  waves  the  automatic  picks  generally  agree  with  the  manual  picks;  at  greater  distances 
the  automatic  picks  are  generally  late.  We  used  this  information  to  construct  a  distance  weighting 
scheme  giving  full  weight  for  P  waves  from  0°  to  2°  distance  with  a  linearly  decreasing  weight 
out  to  8°.  For  S  waves,  full  weight  was  given  up  to  1°  with  a  linear  decrease  to  zero  weight  at  2°. 
However,  event  locations  were  not  very  sensitive  to  changes  in  the  velocity  model  or  the 
weighting.  Locations  based  on  manual  and  on  automatic  picks  are  also  very  similar.  For  the  16 
manually  picked  events,  the  mean  epicenter  difference  is  less  than  3  km  and  for  a  set  of  50 
2.  1<Ml<4.8  events  in  the  PXR  (Monigle,  et  ah,  2012),  we  found  a  median  difference  of  less  than 
2  km  with  a  maximum  difference  of  less  than  8  km  (see  section  6. 1).  The  largest  discrepancies 
occurred  for  events  prior  to  the  completion  of  the  entire  array. 
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Figure  5.  Differences  between  manually  and  automatically  picked  arrival  times  for  16 
stronger  regional  events.  P-differences  in  blue,  S-differences  in  red,  medians  are  solid  lines 
(P:  cyan,  S:  magenta).  Note  the  small  differences  for  a  deep  event  (with  impulsive  P-onset)  for 
all  observed  distances.  Based  on  these  observations,  we  use  S  arrivals  for  event  location  only 
from  0  to  1.5  degrees.  Similarly,  we  downweight  P  waves  arrivals  from  epicentral  distances 
greater  than  3  degrees. 

The  final  seismic  event  database  of  over  22,500  events  (about  50  events  per  day)  passing  all 
location  quality  tests  includes  over  7,900  events  (Figure  6)  with  at  least  25  P  and  S  arrivals, 
corresponding  to  well-located  events  throughout  the  map  area.  Of  these,  4,500  have  a  local 
magnitude  Ml(P)>2.0.  The  locations  based  on  automatic  arrival  time  picks  utilize  S-arrivals, 
constraining  the  epicenters,  and  correlate  with  geologic  structures  seen  in  satellite  imagery 
suggesting  they  are  quite  accurate;  for  some  events  the  USGS  PDE  locations  show  a  mismatch  of 
up  to  50  km  (see  section  6.1  for  quantitative  analysis  of  PDE  and  automatic  location  bias  in  the 
PXR).  Seismicity  is  distributed  throughout  the  study  area.  Most  shallow  seismicity  north  of  the 
Yarlung-Tsangpo  Suture  (YTS),  i.e.,  north  of  29.5°,  is  in  an  area  not  investigated  by  previous 
experiments.  The  most  active  source  regions  are  north-south  trending  grabens  in  the  Lhasa 
terrane.  We  located  almost  10,000  events  in  the  Payang  Basin  (near  30°-31°N,  84°E)  following 
two  Mw>6  main  shocks  in  July  2004  and  April  2005,  respectively.  The  distance  from  the  Payang 
Basin  to  the  closest  stations  is  more  than  100  km  precluding  consideration  for  GT-level  locations 
(in  addition  to  a  one-sided  station  distribution).  However,  our  locations  for  the  two  main  shocks 
are  consistent  with  the  InSAR-derived  fault  area  (Elliot,  et  ah,  2010).  Vigorous  seismicity 
(1,300+  events)  was  also  observed  from  the  vicinity  of  the  PXR  (~30°N,  88°E).  There,  due  to 
close-by  stations,  differences  between  automatic  and  manual  locations  are  only  a  few  kilometers 
(see  section  6.1).  Comparison  with  InSAR  and  manual  locations  suggests  that  epicenter 
uncertainties  for  most  events  in  Figure  6  are  probably  less  then  10  km. 

Adding  requirements  of  distance  to  closest  station  Dmjn<50  km,  azimuthal  gap  Az<225°,  and 
Ml(P)>2.0  results  in  a  high-quality  subset  of  about  600  seismic  events.  In  addition  to  tighter 
clustering,  event  depths  are  well  resolved:  events  north  of  the  southern  Lhasa  Block  are  generally 
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shallow  (<20km)  while  events  south  of  it  are  either  shallow  or  close  to  the  Moho  (Figure  6)  as 
observed  in  previous  studies  (e.g.,  Langin,  et  ah,  2003;  Monsalve,  et  ah,  2006;  Liang,  et  ah,  2008). 
Events  close  to  the  Moho  propagate  across  the  entire  array  with  mantle  velocities  and  therefore 
their  depth  estimates  are  completely  governed  by  the  Moho  depth.  In  tib3p  the  Moho  is  at  7 1  km 
depth  but  in  reality  it  decreases  south  of  the  Yarlung-Tsangpo  Suture  to  about  50  km  under  the 
Himalayas  (Nabelek,  et  al.,  2009).  The  depths  of  the  deep  events  can  be  simply  shifted  to  the 
appropriate  hypocenter  depth  without  altering  the  travel  time  fit  by  choosing  the 


Figure  6.  Left.  Map  of  7,900+  events  with  at  least  25  P  and  S  picks  shows  high-quality 
locations  with  distinct  clustering.  Colors  distinguish  crustal  (yellow)  from  near-Moho  events 
(red;  z>65km).  Black  triangles  are  June  2004  to  August  2005  stations.  PB  is  Payang  Basin, 
PQX Pumqu-Xianza  Rift,  and  YTS  Yarlung-Tsangpo  Suture.  Red  events  outside  network,  e.g. 
in  the  Qiangtang  terrane,  have  poorly  determined  depth.  The  numerous  shallow  events  north 
of  the  Yarlung-Tsangpo  suture  (>29°N)  have  not  been  covered  by  previous  experiments. 
Figure  6.  Right.  Requiring  Dmi„<50  km,  gap  <225°,  and  ML>2.0  results  in  -600  events  with 
well-constrained  depth.  The  depth  distribution  is  bimodal  as  observed  previously;  deep  events 
are  observed  beneath  the  Tethyan  Himalaya  (mainly  between  ~86-88°E)  and  near  the  Yadong- 
Gulu  Rift. 

4.2.  Polarization  S-Picker 


We  used  the  Antelope  software  for  automatic  detection  and  arrival  time  picking  and  for  event 
association  and  location.  The  arrival  picks  were  obtained  through  STA/LTA  analysis  on  vertical 
component  seismic  data.  Our  arrival  time  database  does  contain  a  significant  number  of  S-wave 
picks.  However,  due  to  picking  from  vertical  components,  the  quality  of  the  S  picks  overall  is 
relatively  low  and  degrades  rapidly  with  epicenter  distance  (Figure  5)  forcing  us  to  use  S  arrivals 
for  location  only  to  2°  epicenter  distance.  While  a  few  accurate  S  picks  from  close-by  stations 
generally  suffice  to  constrain  seismic  event  locations,  additional  accurate  S  picks  are  important 
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for  3D  S-wave  velocity  inversions.  Combined  P-  and  S-wave  tomography  provides  critical 
infonnation  regarding  composition  and  physical  state  of  the  crust  and  upper  mantle. 

We  implemented  a  new  automatic,  polarization-based,  S-wave  picker  to  increase  the  number 
of  S  picks  and  to  obtain  accurate  picks  at  distances  larger  than  l°-2°.  The  S  picker  uses  signal 
attributes  from  three-component  seismograms  to  determine  the  S-wave  onset  from  shear  wave 
energy  and  polarization  analysis.  The  picker  has  been  adapted  from  an  implementation  by  Diehl, 
et  al.,  2009  with  polarization  analysis  following  Cichowicz,  1993.  First  we  rotate  the  ZNE  data 
into  the  LQT  ray-coordinate  system  to  separate  P-wave  (L-component)  from  S-wave  energy;  the 
T-component  includes  SH  and  the  Q-component  SV  energy,  respectively.  For  rotation,  we  used 
the  ‘known’  event  location  and  velocity  model,  since  the  events  had  already  been  identified  and 
located  with  picks,  primarily  P  phases,  from  the  Antelope  STA/LTA  algorithm.  The  traces  are 
then  filtered  using  a  Wood- Anderson  and  a  0.5  Hz  high  pass  filter.  The  Wood-Anderson  filter 
integrates  the  velocity  signal  to  displacement,  essentially  acting  as  a  low-pass  filter. 

An  S  phase  is  picked  when  rectilinearity,  directivity  relative  to  the  incoming  P  wave,  as  well 
as  ratio  of  transverse  to  overall  energy  and  the  transverse  amplitude  reach  a  data-dependent 
threshold  within  a  short,  user-defined  time  window.  Rectilinearity  R(t)  is  defined  as  the  degree  of 
linear  polarization  of  the  incoming  particle  motion  (e.g.,  Cichowicz,  1993;  Diehl,  et  al.,  2009); 
R(t)  is  one  for  body  waves.  Directivity  D(t),  which  represents  the  angle  between  the  longitudinal 
(L)  direction  and  the  eigenvector  associated  with  the  maximum  eigenvalue  (see,  e.g.,  Cichowicz, 
1993),  is  zero  for  P  and  one  for  S  waves.  The  energy  ratio  H(t)  is  the  S-wave  energy  (sum  of 
squared  amplitudes)  relative  to  total  energy  and  is  maximum  for  S  arrivals.  The  weight  function 
W(t)  is  simply  the  square  of  the  ratio  of  the  maximum  amplitude  on  either  the  T  or  Q  component 
within  a  short  picking  interval  relative  to  the  overall  maximum  amplitude  on  either  component 
within  the  S-wave  picking  window  (see  below);  we  use  the  squared  amplitude  to  accentuate  the 
S  wave  arrival  (cf.  Diehl,  et  al.  2009).  W(t)  reaches  its  maximum  Wmax  when  the  picking 
interval  includes  the  overall  maximum  amplitude.  We  combine  the  polarization,  R(t)  and  D(t), 
and  amplitude,  H(t)  and  W(t),  attributes  in  a  characteristic  function  CF(t): 

CF(t)  =  R2(t)  *  D2  (t)  *  H2  (t)  *  W(t)  (!) 

CF(t)  is  near  zero  for  incoming  P  wave  energy  and  near  one  for  incoming  S  wave  energy. 

The  procedure  is  illustrated  in  Figure  7  and  variables  and  their  values  are  explained  in  Table  1. 
The  three  traces  in  Figure  7a  show  the  LQT-rotated  seismograms  including  the  P-phase  pick  and 
the  predicted  (gray)  S-phase  arrival  Spre.  The  predicted  arrival  time  is  for  the  initial  ‘known’ 
event  location  and  the  tib3p  velocity  model  (Figure  4).  Figure  7b  shows  the  characteristic 
function  CF(t)  and  the  time  window  from  SW1  to  SW2  within  which  the  S  pick  has  to  occur;  the 
window  is  centered  on  the  predicted  S-arrival  time  and  its  length  depends  on  epicenter  distance. 
The  second  trace  shows  the  weight  function  W(t). 

A  coarse  pick  at  time  WS  (orange  vertical  line  in  Figure  7b)  is  obtained  when  the  weight 
function  W(t)  (blue  line)  exceeds  a  threshold  trl  (thin  red  line)  detennined  from  the  mean  and 
standard  deviation  of  the  weight  function  W(t)  in  the  window  W1  to  W2.  The  window  end,  W2, 
coincides  with  the  earliest  possible  S-arrival  time,  SW1,  and  the  window  start  W1  is  at  the  mid¬ 
point  between  the  P  arrival  and  W2.  This  window  is  chosen  to  exclude  direct  P  energy,  to 
precede  the  S  arrival,  and  to  be  long  enough  to  characterize  the  T  and  Q  component 
‘background’  level.  The  threshold  trl  is  set  sufficiently  higher  than  the  measured  noise  level  so 
that  the  coarse  pick  WS  clearly  identifies  the  S  energy.  (In  this  example,  the  ‘background’  level 
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Figure  7.  Step-by-step  process  for  S  picking  by  polarization  method  for  station  H1400, 108 
km  from  the  3  July  2005  Ml(P)=3.0  earthquake  at  31.55°N  and  85.62°E.  For  details  of  the 
picking  process  refer  to  text.  A)  Rotated  L,  Q,  T  components  with  predicted  S  arrival  Spre.  B) 
Characteristic  function  CF(t)  and  weighting  function  W(t)  built  from  the  seismic  traces.  SW1 
and  SW2  define  picking  window  around  Spre.  Threshold  trl  is  calculated  between  W1  and 
W2;  WS  marks  where  W(t)  crossed  trl.  Wmax  is  the  maximum  amplitude  ofW(t).  C)  WS  is 
projected  on  CF(t).  SWla  is  the  first  minimum  of  CF(t)  prior  to  WS  and  the  threshold  tr2  (red 
line)  is  built  between  SW1  and  SWla.  Zoomed-in  CF(t)  trace  to  right.  D)  Wmax  is  projected 
on  CF(t).  The  final  S  pick  is  the  first  non-local  minimum  below  the  treshold  tr2  on  CF(t)  prior 
to  Wmax.  Zoomed-in  CF(t)  trace  to  right.  E)  Final  S pick  (red)  on  seismic  traces,  zoomed-in 
on  right. 
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of  W(t)  is  very  low  and  therefore  trl  is  also  appears  quite  low).  We  expect  the  actual  S  arrival  to 
closely  precede  WS.  If  WS  cannot  be  found  for  the  given  trl,  no  pick  is  attempted,  thus 
eliminating  data  with  low  signal-to-noise  ratio. 

For  the  final  pick,  we  use  the  characteristic  function  CF(t),  which  incorporates  S  wave  signal 
strength  attributes,  H(t)  and  W(t),  as  well  as  polarization  attributes,  R(t),  D(t)  and  H(t).  For  a  pick 
(Figure  7c),  CF(t)  (blue)  has  to  exceed  a  threshold  tr2  (red  horizontal  line),  defined  from  mean 
and  standard  deviation  of  CF(t)  between  times  SW1  and  SWla,  which  is  the  first  minimum  of 
CF(t)  prior  to  WS.  To  avoid  spurious  maxima  in  CF(t)  prior  to  the  S  arrival,  the  algorithm  works 
backwards  in  time  along  CF(t)  starting  at  Wmax.  A  pick  is  declared  when  CF(t)  reaches  a 
mimimum  after  it  has  dropped  below  the  threshold  tr2,  i.e.,  before  CF(t)  reaches  the  threshold, 
and  stays  below  the  threshold.  The  example  shown  in  Figure  8  shows  a  case  where  CF(t)  had 
reached  a  local  minimum.  To  avoid  local  minima,  the  algorithm,  for  several  samples,  continues 
its  search  backwards  in  time  to  be  able  to  ‘climb’  over  preceding  maxima  and  to  find  the  instance 
where  CF(t)  reaches  a  minimum  after  dropping  below  the  threshold  tr2  and  stays  below  it. 

Figure  8  shows  additional  examples  of  S  picks  as  a  function  of  epicenter  distance,  from  65  to 
255  km,  for  the  upper-crustal  event  of  Figure  7.  All  arrivals  are  crustal  Sg  phases  considering  the 
70  km  thick  crust  of  the  Lhasa  block.  The  picks  are  very  good.  For  this  event,  we  obtained  34  S- 
phase  picks  with  the  polarization  S  picker  compared  with  5  S  picks  from  Antelope  STA/LTA 
analysis.  Figure  9  compares  polarization  picks  to  STA/LTA  based  picks  for  different  distances 
for  another  crustal  seismic  event.  For  the  set  of  common  stations,  the  STA/LTA  picks  are  0.7  s 
later  than  our  new  S  picks.  It  is  no  surprise  that  STA/LTA  based  picks  are  late  considering  the 
higher  noise  level  in  the  P  coda  that  precedes  the  actual  S  arrival;  however  is  some  cases  the 
STA/LTA  picks  can  also  be  early. 

We  have  adjusted  the  picking  parameters  and  thresholds  multiple  times  to  optimize  automatic 
picking  results  for  our  data  based  on  visual  inspection  of  the  algorithm’s  perfonnance  and  on 
comparison  with  manual  picks  from  events  distributed  throughout  south-central  Tibet  (see 
below).  No  pick  is  declared  when  the  thresholds  for  W(t)  or  CF(t)  are  large  or  when  CF(t) 
oscillates  around  its  threshold. 
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Figure  8.  Examples  of  S-phase  picks  with  increasing  distance,  from  65  to  255  km,  for  same 
event  as  in  Figure  7.  Red  are  automatic  S  picks  using  the  new  picker. 


Figure  9.  Polarization  (red)  vs.  STA/LTA  (green)  S  picks  for  the  15  July  2004  ML(P)=2.9 
earthquake  (30.66°N,  83.83°E).  On  average  STA/LTA  S picks  are  0.7 s  late,  although  the 
STA/LTA  picks  can  also  occur  before  the  polarization  S pick  (and  the  actual  arrival).  Because 
the  STA/LTA  picks  were  made  with  P-waves  as  the  primary  target  (vertical  component  and 
higher  (2-10  Hz)  frequency  band),  their  substatially  poorer  performance  is  not  surprising. 
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The  polarization  S  picker  led  to  a  significant  increase  in  the  number  of  successful  S  picks  and 
improved  the  pick  precision  compared  to  STA/LTA  method..  For  example,  for  the  PXR  dataset 
of  50  high-quality  events  (Monigle,  et  al.,  2012),  the  STA/LTA  picker  had  found  478  S  arrivals, 
while  the  polarization  picker  found  1147.  Even  more  important,  at  shorter  distances,  where  pick 
quality  is  high,  the  number  of  S  picks  almost  doubled  from  1 15  to  206.  The  significant  difference 
at  short  distances  (<100  km)  is  critical  for  accurate  hypocenter  locations;  there  the  polarization 
picker  detected  many  more  S  arrivals  and,  on  average,  the  polarization  and  manual  arrival  times 
are  identical. 

Next,  we  applied  the  polarization  S  picker  to  5,790  Ml(P)>2.0  events  within  2°  from  the 
closest  station  and  with  at  least  15  picks  (P  and  S  combined)  in  our  original  database.  These 
events  already  have  well-constrained  epicenter  locations  with  well-recorded  signals  across  a 
significant  part  of  the  Hi-CLIMB  array.  We  do  not  expect  their  locations  to  change  much  when 
changing  and  adding  S  picks.  For  the  5352  shallow  events  (z<30  km),  Sg  can  be  picked  reliably 
to  the  Sg/Sn  crossover  distance  of  about  3°  while  Sn  arrivals  are  either  absent  or  too  weak  to  be 
picked  reliably.  In  contrast,  the  222  deep  events  (z>60  km)  beneath  the  southern  Tibetan  Plateau 
and  the  High  Himalayas  produce  clear  S  arrivals  that  can  be  picked  to  6°-7°  distance,  crucial  to 
constrain  the  uppennost  mantle  S-wave  velocity  structure  beneath  southern  Tibet.  We  obtained 
about  93,000  S  picks  almost  doubling  the  number  of  picks  compared  to  STA/LTA-based  picks 
and,  at  short  distances  of  A<1°,  more  than  tripled  the  number  to  about  6,500  S  picks. 

Reduced  travel  time  plots  for  shallow  and  deep  events  (Figure  10  and  1 1)  show  a  decrease  in 
outliers  and  false  detections  for  the  polarization  S  picks.  For  shallow  events,  S  arrivals  are 
generally  consistent  with  an  average  crustal  velocity  of  3.5  km/s;  however,  as  the  waveform 
examples  show,  the  arrival  time  spread  is  real  and  likely  due  to  velocity  variations  in  the  shallow 
crust.  For  deep  events,  Sn  arrivals  are  present  out  to  about  7°  and  show  a  mantle  S  velocity  of 
about  4.7  km/s.  Large  Sn  arrival  errors  present  beyond  4  degrees  from  the  STA/LTA  analysis 
were  eliminated  through  the  polarization  picking  process. 
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Figure  10.  S-wave  travel  times  vs.  distance  for  small  (2.0<ML(P)<2.7)  seismic  events  shown 
with  3.5  km/s  reduction  velocity.  Polarization  picks  (red)  on  top  ofSTA/LTA  picks  (blue). 
Note  increased  number  of picks  and  reduction  of  outliers  and  false  detections  with 
polarization  picker.  Seismograms  show  spread  in  arrival  times  is  real  indicative  3D  S-wave 
velocity  structure. 


Figure  11.  S-wave  travel  times  vs.  distance  for  deep  (z>60  km)  and  larger  (ML(P)>2.7) 
seismic  events  shown  with  4.7  km/s  reduction  velocity.  Polarization  picks  (red)  on  top  of 
STA/LTA  picks  (blue).  Mantle  S  can  be  picked  consistently  to  6°-7°  distance  with  the 
polarization  picker.  Again  we  obeserve  reduction  of  outliers  and  false  detections  with 
polarization  picker. 

We  relocated  the  5790  seismic  events  with  the  new  polarization  S  picks.  As  expected  for 
events  with  many  phase  picks,  locations  remained  largely  unaffected  (average  horizontal  shift  is 
1 . 1  km)  except  for  events  at  the  periphery  of  the  network;  these  moved  slightly  towards  the 
network.  Depth,  on  average,  changed  by  less  than  one  km  even  though  more  S  picks  were 
available  per  event.  This  is  likely  due  to  a  combination  of  station  coverage  and  number  of  pre¬ 
existing  P  arrivals,  which  yielded  accurate  initial  hypocenters.  In  summary,  out  of  5790  events 
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only  156  changed  depth  by  ±5km  or  more.  Horizontally,  only  85  events  changed  location  by  0. 1 
degrees  or  more.  For  shallow  events  only  1.5%  (88  out  of  5352)  changed  depth  by  ±5km  or 
more.  For  events  at  60km  depth  or  deeper  14.4%  (32  out  of  222)  changed  depth  by  ±5km  or 
more. 

To  further  evaluate  picking  quality  we  compare  relocation  residuals  for  polarization  and 
STA/LTA  arrivals  (Figure  12).  For  shallow  events,  average  residuals  were  -0.07  s  for 
polarization  picks  and  -0.95  s  for  STA/LTA  picks.  For  mid-range  distances,  l°-2.5°,  the  average 
STA/LTA  arrival  residual  is  negative  (arrivals  are  consistently  picked  late).  The  polarization 
picks,  however,  are  centered  on  a  residual  value  of  zero.  The  most  striking  difference  is  at  close 
distances  (<1°),  where  the  number  of  arrivals  increased  while  the  residual  spread  decreased. 


1  1.5  2 

Distance  (degrees) 


Distance  (degrees) 


Figure  12.  Top.  S  residuals  after  event  relocation  with  polarization  picks  (red)  and 
STA/LTA  picks  (blue)  for  shallow  events.  Solid  and  dashed  lines  are  moving  window 
averages  for  polarization  and  STA/LTA  picks  respectively.  STA/LTA  picks  are  consistently 
negative,  meaning  picks  are  late. 

Figure  12.  Bottom,  zoomed-in  for  distances  less  than  1°. 

4.3.  Ml(P) 

We  used  the  Antelope  module  dbml  to  calculate  local  magnitude  (ML)  from  the  vertical 
component  seismograms  for  all  events  in  the  database  with  epicenter-station  distances  of  600  km 
or  less.  Usually  ML  is  determined  from  the  maximum  amplitudes  on  Wood- Anderson-equivalent 
horizontal  seismograms  with  the  maximum  amplitude  corresponding  to  an  S-phase.  We  thus 
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originally  defined  the  time  window  to  be  1.5  times  the  S-P-time  for  a  given  station-event 
distance  to  obtain  the  amplitudes.  However,  during  times  of  vigorous  seismicity,  we  often 
encountered  contamination  by  arrivals  from  closely  following,  larger  events  (Figure  13).  Another 
complication  arises  from  using  vertical  component  data  that  underestimates  S  arrival  amplitude, 
which  led  to  systematically  underestimated  magnitudes. 

To  avoid  both  problems,  we  decided  to  estimate  event  size  from  a  short  time  window 
containing  only  P  waves  (0.25  S-P-time,  Figure  13).  This  approach  eliminated  contamination  and 
provided  a  much  cleaner  magnitude  dataset.  Using  the  average  theoretically  expected  S-to-P 

•5 

amplitude  ratio  of  (vp/vs)  ~  5  for  earthquakes  with  a  double-couple  mechanism  provides  a 
physics-based  approach  to  convert  to  S-based  magnitude  by  simply  adding  a  constant  value  of 
0.7.  We  call  this  (corrected)  magnitude  Ml(P)  and  use  it  consistently  throughout  this  report. 

On  average,  ML(P)  is  consistent  with  moment  magnitude  Mw  from  regional  moment  tensors 
(Baur,  2007)  and  matches  teleseismic  Global  CMT  Mw  and  USGS  NEIC  body  wave  magnitude 
mb  (Figure  14). 

Our  event  Ml(P)  are  medians  based  on  at  least  3  observations  within  600  km  distance  and  a 
minimum  signal-to-noise  ratio  of  3  after  Wood- Anderson  filtering  resulting  in  over  4,000  events 
with  Ml(P).  The  frequency-magnitude  plot  suggests  completeness  for  ML(P)>3.3  within  600  km 
and  Ml(P)>2.0  within  150  km  of  the  array  (Figure  15). 


MI190  DHZ 

Ml  700  BH/ 
MI 750  BHZ 
Ml  7  30  8  HZ 
Ml  TOO  8 HZ 
HIZ40  8 HZ 
HUGO  8 HZ 
HUSO  BH7 
HIZ70  8 HZ 
HUSO  BHZ 
HI 780  BHZ 
Ml  170  BHZ 
HI110  BHZ 
HI  310  BHZ 
Ml  471  BHZ 
H0810BHZ 
HOOOO  BHZ 
HI 370  8 HZ 
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Figure  13.  Wood-Anderson  equivalent  waveforms  for  90  seconds  of  data  from  April  8, 
2005,  about  8  hours  after  an  M=6+  earthquake  in  the  Payang  Basin.  A  small  event  (first  set 
of  arrivals)  is  followed  by  a  much  stronger  event  (second  set  of  P  arrivals)  within  1.5  times  its 
S-P-time  (purple  area).  The  MLfor  the  first  event  based  on  the  long  time  window  is  about  2 
magnitude  units  too  large.  Another  problem,  apparent  for  the  second  event,  is  that  vertical 
seismograms  underestimate  the  amplitude  of  S-phases,  which  are  predominantly  recorded  on 
horizontal  components;  in  this  case  the  S-to-P  amplitude  ratio  is  about  2  compared  to  the 
expected  value  of  ~5.  We  decided  to  determine  magnitudes  from  a  short  window  containing 
only  the  P  arrival  (0.25  S-P-time  shown  in  brown). 
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Figure  14.  Magnitude  comparison.  ML  is  local  magnitude  from  P  wave  amplitudes. 
MW(RMT),  MW(CMT),  nth(NEIC)  are  moment  magnitude  from  regional  moment  tensor 
inversion,  Global  CMT  analysis,  and  USGS  NEIC  body  wave  magnitude,  respectively. 
Shown  are  common  events.  The  average  underestimation  of  uncorrected  ML  of  0.5,  0. 7,  and 
0. 7,  respectively,  is  consistent  with  the  expected  theoretical  bias  due  the  use  of  P  instead  of 
S  waves.  The  scatter  around  the  average  (short  dashes,  long  dashes  correspond  to  one-to- 
one  correspondence)  is  similar  to  scatter  observed  when  regressing  other  magnitudes 
relative  to  each  other  (e.g.,  Braunmiller,  et  al.,  2005). 
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Figure  15.  Left.  ML(P)  versus  NEIC  mb,  Global  CMT  and  regional  moment  tensor  (Baur, 
2007)  Mw. 

Figure  15.  Right.  Frequency-magnitude  relation  for  4,000+  events  within  600  km  of  the 
array  and  at  least  3  magnitude  observations  with  minimum  signal-to-noise  ratio  of  3  per 
event. 


5.  SOURCE  PARAMETERS  FROM  REGIONAL  MOMENT  TENSORS 

Regional  moment  tensor  analyses  was  performed  (Burtin,  2005;  Burtin,  et  al.,  2005;  Baur, 
2007;  Monigle,  et  al.,  2012)  for  143  selected  larger  earthquakes  in  the  Himalayan-Tibetan  region 
using  the  full  wavefonn  inversion  code  developed  by  Nabelek  and  Xia,  1995.  Results  are 
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summarized  in  Figure  16  and  listed  in  Table  1.  Our  efforts  substantially  expand  the  region’s 
moment  tensor  catalog;  our  database  includes  most  events  analyzed  by  de  la  Torre,  et  al,  2007. 

Burtin,  2005  tested  analysis  feasibility  with  permanent  far-regional  stations.  Deriving 
station  specific  velocity  models,  he  analyzed  27  Mvv>4.5  events  that  occurred  from  1988  to  2001 
within  1000  km  of  the  closest  station  LSA  (Lhasa,  Figure  1)  and  along  the  southern  boundary  of 
the  Tibetan  Plateau.  Long  paths  require  analysis  at  long  periods  (T>33  s),  which  limits  parameter 
resolution  and  sets  a  high  analysis  threshold.  For  seismic  events  during  1991-1992,  he  also  used 
data  from  the  temporary  Passcal  91-92  deployment  (e.g.,  Zhu  and  Flelmberger,  1996). 

79  80  81  82  83  84  85  86  87  88  89  90  91  92  93  94  95  96  97  98 


79  80  81  82  83  84  85  86  87  88  89  90  91  92  93  94  95  96  97  98 


0  2000  4000  6000 


Elevation 

Figure  16.  Regional  moment  tensor  solutions.  (From  Monigle,  et  al.,  2012  and  after  Burtin, 
2005  and  Baur,  2007.)  We  analyzed  143  seismic  events  in  Nepal  and  southern  Tibet  using 
high-quality  broadband  waveforms.  The  Hi-CLIMB  set  of  115  solutions  contains  events  as 
small  as  Mw=3.5.  In  southern  Tibet  and  Nepal  our  efforts  substantially  expand  the  moment- 
tensor  catalog. 


The  Hi-CLIMB  data  allow  analysis  of  smaller  Mw>3.5  events  with  near-regional 
seismograms  at  periods  of  T>15-20  s  (Baur,  2007;  Monigle,  et  al.,  2012).  We  obtained  115 
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moment  tensors.  High  quality  and  redundancy  of  the  data  allow  parameter  resolution  tests  for 
size,  depth,  and  moment  tensor.  For  the  PXR  events  (near  30°N,  88°E),  where  we  have  good 
event  location  control  (Monigle,  et  ah,  2012),  a  significant  number  of  events  with  similar-looking 
waveforms  spanning  a  wide  size  range  and  good  recording  geometry  allowed  testing  of 
parameter  resolution  and  evaluating  of  velocity-depth  models  for  seismogram  calculation.  The 
tests  suggest  a  depth  resolution  for  shallow  crustal  events  on  the  order  of  ±5  km,  an  estimated 
Mw  uncertainty  of  ±0. 1-0.2  units,  and  uncertainties  for  strike,  dip,  and  rake  on  the  order  of  5°- 
10°  with  currently  used  velocity  models  and  event  locations. 

Table  1. 143  regional  moment  tensor  solutions  in  Nepal  and  southern  Tibet.  S/D/R  are 
strike,  dip  and  rake  of  one  nodal  plane;  CD  is  centroid  depth;  Mw  is  moment  magnitude,  DC  is 
double-couple  percentage  of  deviatoric  moment  tensor;  ST  is  number  of  stations  used  for 
inversion;  Aut  designates  author,  Bu:  Burtin,  2005,  Ba:  Baur,  2007,  Mo:  Monigle,  et  al., 

2012. 


Event  ID 
YrMoDa  HrMi 

Lat 

(°) 

Lon 

(°) 

S/D/R 

(77°) 

CD 

(km) 

Mw 

DC 

(%) 

ST 

(#) 

Aut 

880820  2309 

26.755 

86.616 

317/89/-1 1 1 

52 

6.7 

100 

2 

Bu 

881029  0910 

27.871 

85.648 

284/53/87 

21 

5.1 

53 

2 

Bu 

911209  0102 

29.543 

81.632 

114/53/83 

17 

5.0 

64 

9 

flu 

911221  1952 

27.795 

87.955 

212/88/13 

66 

4.9 

57 

4 

an> 

920206  0335 

29.61 

95.521 

204/63/87 

11 

5.2 

47 

6 

mm 

920404  1743 

28.147 

87.979 

145/66/-153 

55 

5.0 

41 

7 

wm 

920602  2207 

28.984 

81.913 

346/79/-151 

53 

5.2 

61 

10 

11 

920730  0824 

29.584 

90.163 

1 74/60/- 1 17 

9 

6.1 

90 

2 

am- 

930320  1451 

29.084 

87.333 

359/53/-95 

13 

6.3 

67 

3 

■a 

930331  1344 

29.091 

87.349 

336/6 1/-124 

16 

5.1 

85 

3 

Bu 

930524  0502 

28.835 

96.082 

87/72/90 

22 

4.9 

55 

7 

mm 

950217  0244 

27.635 

92.371 

238/76/18 

25 

5.3 

100 

6 

B3K 

960426  1630 

27.825 

87.821 

84/79/- 112 

75 

4.7 

93 

1 

wm 

960609  2325 

28.325 

92.201 

210/65/-85 

80 

5.3 

79 

4 

mm 

960925  1741 

27.433 

88.552 

140/70/105 

25 

4.8 

70 

5 

mm 

961120  2327 

28.853 

96.021 

321/55/83 

33 

4.7 

14 

2 

wm 

970105  0847 

29.845 

80.532 

108/71/79 

17 

5.4 

64 

6 

an 

970718  1939 

26.811 

91.793 

7/66/139 

27 

4.5 

93 

3 

Bu 

971030  0202 

29.552 

89.698 

143/87/- 142 

94 

5.3 

82 

3 

Bu 

971103  0229 

29.078 

85.383 

1 62/48/- 1 10 

11 

96 

7 

m 

980708  0344 

27.325 

73/85/101 

71 

6 

mm 

980818  0410 

27.55 

121/87/85 

21 

95 

4 

1:J 

980926  1827 

27.77 

60/58/95 

16 

86 

4 

an 

981126  1014 

27.753 

19/5  8/- 106 

57 

62 

1 

mm 

000102  1023 

27.559 

92.498 

74/75/77 

15 

4.8 

57 

3 

Bu 

000125  1643 

27.663 

92.631 

69/67/71 

13 

4.8 

67 

2 

Bu 

010927  2240 

26.817 

87.593 

50/54/-53 

27 

3.8 

75 

5 

Ba 

011106  1409 

27.268 

91.919 

79/68/87 

15 

4.5 

72 

12 

Ba 
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011107  0458 

29.99 

85.6* 

349/54/-101 

22 

4.1 

84 

9 

Ba 

011127  0731 

29.606 

81.752 

109/77/85 

15 

5.3 

64 

4 

Bu 

011202  2241 

27.15 

88.171 

117/56/124 

15 

4.8 

89 

13 

Ba 

011218  0222 

29.612 

81.78 

135/68/96 

21 

4.0 

72 

6 

Ba 

020307  1550 

29.56 

84.83* 

137/89/-121 

77 

4.2 

97 

9 

Ba 

020323  0456 

29.863 

87.928 

189/57/-65 

7 

4.5 

72 

11 

Ba 

020402  1957 

29.287 

87.009 

221/80/-43 

18 

3.5 

47 

9 

Ba 

020508  1756 

28.659 

86.584 

232/73/30 

77 

3.8 

97 

6 

Ba 

020509  0109 

29.826 

87.837 

182/53/-67 

6 

4.1 

78 

13 

Ba 

020604  1436 

30.595 

81.44 

326/82/- 1 17 

15 

5.6 

99 

6 

Ba 

020604  1727 

30.425 

81.398 

324/75M6 

27 

4.3 

48 

10 

Ba 

020606  0135 

30.457 

81.251 

306/87/-120 

18 

4.2 

88 

3 

Ba 

020620  0540 

26.25 

89* 

148/85/162 

30 

4.6 

100 

13 

Ba 

020702  1710 

27.11 

84.78* 

173/74/- 148 

70 

3.9 

66 

6 

Ba 

020709  0308 

29.889 

87.995 

325/61/-1 12 

6 

4.4 

94 

10 

Ba 

020716  1839 

27.774 

87.468 

192/58/-69 

65 

3.8 

77 

7 

Ba 

020718  2324 

28.2 

87.927* 

356/87/-37 

78 

3.7 

48 

6 

Ba 

020806  2330 

29.912 

87.87 

332/47/- 108 

6 

4.4 

86 

13 

Ba 

020822  0450 

29.818 

88.087 

32 1/7 1/- 1 07 

9 

4.3 

93 

16 

Ba 

020831  1240 

29.92 

88.12* 

328/70/- 106 

5 

4.6 

91 

23 

Ba 

020927  1714 

33.364 

93.567 

254/79/-25 

39 

5.1 

85 

10 

Ba 

021019  0724 

35.654 

93.129 

2 10/89/- 150 

36 

4.9 

89 

3 

Ba 

021026  2028 

35.141 

96.097 

183/72/141 

21 

5.4 

81 

8 

Ba 

021104  1059 

28.855 

82.208 

262/59/-122 

24 

3.5 

99 

8 

Ba 

021116  0852 

29.667 

90.409 

225/64/-59 

16 

4.4 

83 

15 

Ba 

021129  1649 

29.676 

90.395 

235/64/-54 

16 

4.7 

97 

25 

Ba 

030116  1136 

29.8 

88.006 

339/48/- 1 13 

3 

4.8 

61 

11 

Ba 

030116  2215 

32.615 

89.948 

227/89/-29 

25 

5.0 

95 

10 

Ba 

030118  1031 

28.622 

81.95 

120/58/117 

33 

3.8 

76 

7 

Ba 

030131  0604 

29.746 

90.515 

240/69/-49 

10 

4.0 

74 

9 

Ba 

030211  1036 

32.512 

93.787 

98/86/- 164 

27 

4.8 

50 

2 

Ba 

030226  1958 

28.347 

86.449 

180/72/- 160 

75 

3.9 

93 

8 

Ba 

030325  1851 

27.264 

89.331 

229/87/23 

32 

5.5 

100 

13 

Ba 

030331  0532 

26.528 

91.827 

307/74/45 

65 

3.9 

58 

4 

Ba 

030520  1834 

32.681 

93.094 

226/82/-37 

21 

5.1 

70 

17 

Ba 

030524  1127 

32.588 

92.425 

214/56/-79 

12 

4.5 

82 

2 

Ba 

030524  1932 

32.609 

92.337 

334/86/174 

36 

4.9 

75 

9 

Ba 

030529  1418 

35.705 

80.653 

110/50/79 

21 

4.8 

39 

4 

Ba 

030707  0655 

34.606 

89.475 

327/80/- 168 

21 

5.8 

60 

14 

Ba 

030708  1230 

26.905 

87.546 

222/79/-53 

12 

3.7 

15 

12 

Ba 

030728  0245 

28.698 

82.045 

108/82/79 

15 

3.7 

21 

10 

Ba 

030818  0903 

29.573 

95.605 

250/90/-45 

24 

5.7 

46 

16 

Ba 

030929  1340 

27.358 

87.738 

69/71/45 

12 

4.0 

66 

15 

Ba 

031028  0231 

32.446 

92.173 

217/76/-70 

7 

4.8 

39 

16 

Ba 

031122  0507 

28.277 

83.767 

265/88/-94 

16 

3.8 

87 

22 

Ba 

19 
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031122  2331 

28.202 

83.767 

71/82/81 

21 

3.5 

84 

29 

Ba 

031123  1915 

28.297 

83.773 

270/77/-87 

14 

3.7 

89 

17 

Ba 

031210  1717 

30.469 

83.278 

307/82/- 108 

30 

4.3 

94 

11 

Ba 

031212  1339 

37.626 

94.588 

59/81/24 

25 

4.8 

98 

18 

Ba 

040103  1314 

27.739 

85.931 

109/75/87 

14 

4.7 

96 

25 

Ba 

040106  0313 

30.539 

86.042 

18/83/- 16 

22 

4.6 

97 

4 

Ba 

040210  0539 

32.745 

83.301 

333/65/- 1 10 

29 

4.6 

24 

11 

Ba 

040218  0123 

27.361 

87.758 

291/82/160 

24 

3.9 

75 

23 

Ba 

040224  2021 

37.487 

96.828 

144/80/- 166 

17 

4.9 

57 

19 

Ba 

040227  1253 

28.127 

87.66 

108/54/- 133 

92 

4.6 

80 

21 

Ba 

040301  1741 

30.368 

80.496 

241/65/-128 

23 

4.2 

93 

19 

Ba 

040306  1154 

33.291 

91.953 

3 1/56/- 108 

5 

4.7 

67 

5 

Ba 

040307  1329 

31.639 

91.236 

228/80/-54 

11 

5.6 

81 

9 

Ba 

040316  2123 

37.558 

96.668 

151/69/151 

11 

4.9 

35 

5 

Ba 

040327  1847 

33.954 

89.179 

26/51/-60 

8 

5.8 

93 

11 

Ba 

040327  2005 

33.931 

89.292 

353/63/-126 

9 

4.4 

83 

5 

Ba 

040328  2205 

34.132 

89.275 

324/76/- 159 

23 

4.9 

41 

7 

Ba 

040328  2227 

33.957 

89.265 

337/68/- 130 

15 

4.8 

29 

7 

Ba 

040403  0251 

29.846 

81.118 

309/62/127 

74 

4.1 

77 

4 

Ba 

040422  1002 

33.998 

89.218 

26/53/-56 

5 

5.2 

50 

10 

Ba 

040504  0504 

37.506 

96.758 

56/83/-33 

28 

5.4 

89 

6 

Ba 

040510  2327 

37.485 

96.6 

56/59/66 

12 

5.6 

71 

14 

Ba 

040523  0738 

34.079 

89.285 

229/78/13 

16 

5.3 

78 

7 

Ba 

040523  1446 

34.08 

89.32 

233/88/16 

20 

5.0 

65 

9 

Ba 

040605  0847 

29.86 

89.698 

208/62/-70 

9 

3.9 

72 

3 

Ba 

040624  1003 

29.868 

87.905 

14/54/-76 

6 

4.1 

54 

14 

Ba 

040630  1533 

30.017 

88.097 

183/32/-79 

12 

4.2 

83 

54 

Mo 

040701  1914 

30.016 

88.082 

177/30/-92 

12 

3.5 

95 

24 

Mo 

040703  1410 

34.093 

89.349 

2 04/64/- 80 

10 

5.4 

83 

14 

Ba 

040708  2149 

30.032 

88.161 

1 9 1/33/-70 

12 

4.1 

99 

21 

Mo 

040711  2308 

30.694 

83.672 

348/49/-90 

16 

6.2 

47 

17 

Ba 

040712  1438 

30.746 

83.773 

31/80/-39 

25 

4.6 

42 

23 

Ba 

040715  2339 

30.026 

88.173 

183/30/-81 

12 

3.4 

90 

45 

Mo 

040716  2033 

28.232 

83.944 

108/81/92 

10 

3.6 

73 

14 

Ba 

040720  0335 

27.938 

85.793 

94/73/95 

13 

3.6 

89 

12 

Ba 

040721  1544 

30.035 

88.150 

183/29/-79 

12 

3.7 

83 

58 

Mo 

040723  0125 

30.013 

88.155 

179/30/-78 

12 

4.5 

84 

83 

Mo 

040728  2222 

30.714 

83.633 

2/70/-87 

10 

5.1 

96 

9 

Ba 

040804  0209 

25.923 

90.262 

353/89/54 

53 

4.1 

26 

18 

Ba 

040821  0907 

30.037 

88.159 

177/30/-85 

12 

4.5 

94 

81 

Mo 

040824  1005 

32.542 

92.19 

326/83/165 

12 

5.3 

61 

10 

Ba 

040903  1024 

30.033 

88.150 

176/25/-77 

12 

3.5 

51 

23 

Mo 

040909  0555 

29.466 

80.281 

147/55/88 

6 

3.8 

18 

12 

Ba 

040912  0848 

29.502 

81.575 

127/56/121 

21 

3.9 

94 

10 

Ba 

040915  1310 

30.045 

88.159 

180/26/-78 

12 

3.8 

77 

56 

Mo 
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040923  1728 

31.116 

96.308 

94/80/- 16 

31 

4.4 

92 

9 

Ba 

040927  1705 

29.773 

95.515 

311/76/160 

23 

5.0 

86 

13 

Ba 

041001  1803 

30.046 

88.156 

175/22/-94 

15 

3.4 

57 

32 

Mo 

041026  0211 

31.024 

81.154 

187/62/-35 

13 

5.4 

36 

14 

Ba 

041026  1113 

31.754 

92.893 

76/89/- 16 

25 

4.7 

67 

27 

Ba 

041104  1400 

30.056 

88.161 

164/32/-97 

12 

3.7 

79 

45 

Mo 

041110  0421 

27.929 

87.778 

263/69/79 

66 

4.4 

52 

12 

Ba 

041124  2235 

27.329 

90.944 

76/88/92 

10 

3.6 

66 

5 

Ba 

050108  1251 

29.925 

88.338 

96/79/- 179 

15 

3.2 

87 

11 

Mo 

050115  2232 

29.439 

80.965 

263/62/80 

22 

3.7 

84 

13 

Ba 

050116  0843 

29.646 

80.696 

140/71/101 

10 

4.3 

76 

26 

Ba 

050208  0151 

27.82 

86.033 

99/57/50 

20 

3.8 

98 

22 

Ba 

050208  0335 

27.738 

86.092 

108/75/89 

10 

4.1 

44 

30 

Ba 

050208  0413 

27.76 

86.072 

325/59/102 

21 

3.9 

73 

10 

Ba 

050227  1832 

25.371 

91.553 

332/73/69 

19 

3.5 

77 

8 

Ba 

050407  0140 

29.299 

81.451 

300/74/- 107 

23 

3.6 

78 

10 

Ba 

050407  2004 

30.491 

83.662 

349/56/-95 

8 

6.3 

79 

21 

Ba 

050407  2141 

30.588 

83.681 

355/52/-91 

10 

5.3 

98 

28 

Ba 

050408  1045 

29.977 

83.012 

226/726-85 

10 

4.0 

76 

10 

Ba 

050409  0920 

30.494 

83.599 

230/61/59 

34 

4.4 

49 

15 

Ba 

050508  1642 

30.097 

90.261 

3 17/66/- 149 

98 

4.8 

97 

17 

Ba 

050511  0412 

30.671 

85.244 

32 1/8 1/-166 

4 

4.4 

82 

70 

Mo 

050515  1921 

31.613 

92.676 

55/67/52 

35 

4.2 

91 

16 

Ba 

050522  1417 

30.319 

83.822 

25/88/46 

12 

4.0 

95 

10 

Ba 

050601  2006 

28.881 

94.626 

283/88/-96 

17 

5.5 

48 

14 

Ba 

6.  PUBLISHED  RESULTS  AND  ONGOING  RESEARCH 
6.1.  Low  Angle  Normal  Faulting  in  the  PXR 

Seismicity  north  of  the  Yarlung-Tsangpo  Suture  in  the  Lhasa  Terrane  of  southern  Tibet 
occurs  primarily  along  N-S  trending  grabens.  Analysis  concentrated  on  seismicity  along  the 
Pumqu-Xianza  Rift  (PXR)  during  July  2004  to  June  2005  when  close  sites  T0370  and  T0320 
were  operating  (Figure  17).  The  database  contains  about  350  events  with  25  or  more  automatic 
picks  slightly  west  of  the  PXR  compared  to  1 1  events,  ranging  in  magnitude  (mb)  from  4.0  to 
4.8,  listed  in  the  USGS  NEIC  catalog.  We  relocated  the  40  largest  events  (ML(P)>3.2)  using 
manual  P  and  S  picks.  Our  database  also  includes  events  near  T0370  and  T0320  (Figure  17).  We 
selected  10  stronger  events  (2.1<ML(P)<3.4)  to  verify  their  locations;  six  satisfy  GT5-criteria 
(Table  2,  marked  “&”). 

A  detailed  description  of  the  analysis,  including  manual  earthquake  locations,  regional 
moment  tensors,  and  a  discussion  of  the  mechanics  involved  in  30°  low-angle  normal  faulting 
has  been  published  in  Monigle,  et  ah,  2012.  Here  we  report  only  on  earthquake  location  and 
source  aspects  relevant  to  this  project. 
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Automated  locations  with  local  Hi-CLIMB  stations  reveal  a  systematic  west-southwest  bias 
of  routine  NEIC  epicenters  that  can  reach  ~50  km  (Figure  17).  Hi-CLIMB  locations  based  on 
automatic  P  and  S  picks  are  closer  to  the  rift  and  cluster  near  the  central  PQX.  Manual  phase 
picking  for  selected  events  verified  the  high  automatic  location  quality  with  epicenter  differences 
of  only  a  few  kilometers  (Figure  17). 


Figure  17.  Left.  Comparison  of  NEIC  (green  circles)  and  automated  Hi-CLIMB  (blue 
squares)  seismic  event  locations  for  the  PXR  reveals  a  general  SW  NEIC  location  bias  that 
can  reach  more  than  40  km  (average  value  for  11  common  events  is  25  km). 

Figure  17.  Right.  Automatic  (blue)  and  manual  (orange)  locations  for  50  PXR  events  differ 
on  average  by  less  than  2  km  and  differences  do  not  exceed  8  km.  The  main  northern  cluster 
lacks  azimuthal  coverage  to  satisfy  GT5  criteria. 


The  manual  locations  of  the  main  northern  group  cluster  tightly  west  of  the  central  PXR 
(Figure  17  and  18)  and  at  about  10  km  depth.  The  hypocenters  are  constrained  by  phase  picks 
from  close  stations  and  are  consistent  with  ~12  km  depth  (uncertainties  of  about  ±3  km) 
estimated  from  moment  tensor  analysis  (Monigle,  et  ah,  2012).  All  main  cluster  mechanisms  are 
consistent  with  N-S  oriented  normal  faulting  on  a  -30°  west-dipping  fault.  This  cluster 
technically  does  not  satisfy  GT5  criteria  (e.g.,  Bondar,  et  al.,  2004),  even  though  the  closest  site  is 
-20  km  from  the  epicenters,  because  of  the  large  -150°  azimuthal  gap  since  the  events  occurred 
at  the  northern  terminus  of  the  lateral  array. 

The  seismic  events  south  of  T0370  are  within  the  network  and  satisfy  the  azimuthal  and 
close-distance-recording  criteria  for  GT5;  however,  the  events  are  smaller,  ultimately  only  6 
satisfy  the  distance  requirement.  The  largest  event  of  the  cluster  just  south  of  T0370  has  a  strike- 
slip  mechanism  and  its  waveform  modeling  depth  is  consistent  with  its  hypocenter  depth. 
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88° 


Figure  18.  Map  of  the  central  PXR  showing  manual  locations  for  50  events.  Blue  are  events 
west  of  the  rift,  red  near  station  T0370  and  black  south  of  T0320.  White  dots  are  epicenters 
from  automatic  picks  (Carpenter,  et  al.,  2010).  Green  triangles  are  Hi-CLIMB  stations. 
Diamond  is  town  of  Qingdu.  Open  circles  are  NEIC  locations  for  6  events  that  all  relocated 
into  the  blue  cluster.  Fault  plane  solutions  from  moment  tensor  analysis  are  given  by  event 
date  (YYMMDD).  After  Monigle,  et  al.,  2012. 


6.2.  Velocity  and  Attenuation  Tomography 

In  addition  to  the  immediate  deliverables  of  this  project,  our  current  work  is  focused  on 
obtaining  3D  P-  and  S-wave  velocity  structure  from  inversion  of  local  and  near-regional  body 
waves  and  on  inverting  for  the  crust  and  upper  mantle  attenuation  properties  (Q).  Both  projects 
are  perfonned  by  graduate  students  and  planned  to  be  completed  in  the  near  future. 

6.2.1  Velocity  Tomography 

Resolving  3D  P-  and  S-  wave  velocity  and  structure  of  the  south-central  Tibetan  Plateau  is 
the  focus  of  a  tomography  effort  aimed  at  complementing  and  expanding  on  the  existing  receiver 
function  image  of  the  collision  zone  (Nabelek,  et  al.,  2009).  We  are  using  the  tomoDD  double- 
difference  tomography  code  (Zhang  and  Thurber,  2003)  and  have  begun  with  synthetic  data 
inversions  to  evaluate  how  well  our  existing  P  and  S  wave  dataset  is  able  to  resolve  3D  crust  and 
upper  mantle  structure  considering  the  actual  station  and  event  distribution. 

Our  dataset  consists  of  the  P  and  S  picks  from  the  5790  seismic  events  described  in  section 
4.2  for  which  we  have  at  least  15  P  or  S  arrivals,  which  occurred  within  2°  from  the  closest  Hi- 
CLIMB  station,  and  have  a  minimum  magnitude  of  ML(P)>2.0  to  ensure  significant  path  lengths. 
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These  criteria  result  in  more  than  200,000  P  picks  and,  due  to  the  newly  implemented 
polarization  S  picker,  in  more  than  90,000  S  arrivals  providing  dense  ray  coverage  over  the  study 
area  to  invert  for  a  high-resolution  regional  velocity  structure.  We  are  particularly  excited  about 
our  S  arrival  dataset,  which  allows  a  joint  inversion  for  P-  and  S-wave  structure  providing 
important  compositional  constraints. 

6.2.2  Attenuation  Tomography 

Resolving  the  3D  attenuation  structure  (Q)  of  the  south-central  Tibetan  Plateau  is  the  focus  of 
a  second  tomography  project.  Estimating  3D  Q  structure  from  t*  measurements  is  widely  applied 
(e.g.,  Hauksson  and  Shearer,  2006).  Our  initial  efforts  had  focused  on  estimating  the  resolving 
power  of  our  dataset  to  invert  event-station  specific  t*  values  for  3D  Q  structure  considering 
limited  ray  coverage  in  the  lower  part  of  the  ~70-km  thick  crust  and  in  the  uppennost  mantle 
(through  Pn  and  mantle  P)  beneath  Tibet.  Currently,  we  are  working  on  optimizing  unbiased  t* 
estimation  from  windowed  body  wave  phases;  work  is  perfonned  using  synthetic  (spikes  as  well 
as  full  waveform  synthetic  seismograms)  and  real  data.  Once  we  obtain  good  t*  estimates, 
inverting  for  Q  will  be  fast  since  the  tomography  procedures  are  already  in  place. 

7.  GT5  DATASET 

We  applied  the  Bondar,  et  ah,  2004  criteria  to  identify  candidate  events  for  the  seismic  event 
ground  truth  (GT)  database.  54  events  from  the  5790  event  database  described  in  section  4.2, 
which  is  based  on  STA/LTA  P  picks  and  polarization  S  picks,  satisfied  the  GT595%-level  criteria 
and  were  above  an  ML(P)>2.0  threshold.  Arrival  times  for  all  54  events  were  manually  re-picked 
from  three-component  seismograms  and  then  the  events  were  relocated.  We  applied  the  P 
distance  weighting  described  in  section  4.1  to  both  P  and  S  phases.  The  location  results  contain 
5 1  events  recorded  by  10  or  more  stations  within  250  km  distance,  with  an  azimuthal  gap  of  1 10° 
or  less,  at  least  one  station  within  30  km  from  the  epicenter,  and  a  secondary  azimuthal  gap  of 
160°  or  less.  After  manual  relocation,  a  small  shift  in  the  epicenter  resulted  in  two  events  not 
meeting  the  azimuthal  gap  criterion  and  one  event  was  found  to  be  much  smaller  than  Ml(P)=2.0 
and  could  not  be  picked  reliably.  Of  the  5 1  events,  41  are  at  GT5g5%-level  with  records  extending 
as  P/Pn  beyond  the  250  km  distance.  The  remaining  10  events  are  generally  smaller;  9  of  them 
satisfy  GT595%-level  criteria  but  lack  distant  detections  and  one  event  lacks  distant  arrivals  and 
does  not  satisfy  the  secondary  azimuthal  requirement.  All  51  events  are  shown  in  Figure  19. 

Table  2  includes  all  41  GT595o/o  (candidate)  events. 

For  event  #19,  moment  tensor  inversion  resulted  in  a  double-couple  strike-slip  source 
mechanism  at  15  km  depth  with  roughly  N-S  and  E-W  trending  nodal  planes  (Monigle,  et  ah, 
2012).  Parameter  uncertainties  (-20°  for  strike,  dip,  and  rake  and  ±10  km  for  depth)  are  large  due 
to  small  event  size  (Mw=3.2,  ML(P)=3.1).  The  only  event  listed  in  the  USGS  NEIC  and  the  ISC 
catalogs  is  event  #32.  The  NEIC  and  ISC  locations  are  about  13  km  and  over  30  km, 
respectively,  south-southwest  of  our  location.  Manual  locations  with  Hi-CLIMB  data  (ours  and 
Griffin,  et  ah,  2011)  are  within  a  few  km  of  the  automated  location.  Moment  tensor  analysis 
reveals  a  well-constrained  double-couple  strike-slip  mechanism  at  4  km  depth  with  an  E-W 
trending  ci3-axis.  Magnitude  estimates  for  the  event  are  similar  (Ml(P)=4.7,  Mw=4.4  vs. 
mb(NEIC)=4.6  and  mb(ISC)=4.3). 
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Figure  19.  Map  of  41  GT595o/o  candidate  events  (circles).  Diamonds  are  10  other  candidates 
(all  lack  P/P„  beyond  250  km  distance;  all  except  one  satisfy  the  secondary  gap  criterion). 
Shown  are  potentially  monitoring  relevant  events  with  ML(P)  >2.0  that  can  be  detected  several 
hundred  km  from  a  seismic  event  source.  Color  as  in  Figure  6. 
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Table  2.  Manual  locations  of  41  GT5<)5%  (candidate)  events.  are PXR  events.  Moment 


tensor  analysis  was  performed  for  events  #19  and  #32. 


Event 

(#) 

Date 

(M/D/Y) 

Time 

(UTC) 

Latitude 

(°N) 

Longitude 

(°E) 

Depth 

(km) 

ml(P) 

1 

06/17/04 

07:00:23.7 

31.607 

85.066 

6.3 

2.2 

2 

06/26/04 

01:26:36.7 

31.590 

85.026 

0.7 

2.0 

3 

07/02/04 

00:45:10.7 

31.608 

85.043 

0.7 

2.0 

4 

07/02/04 

01:00:09.4 

31.611 

85.057 

0.9 

2.3 

5 

07/22/04 

17:34:29.7 

28.403 

88.129 

72.0 

2.1 

6 

08/02/04 

21:48:29.6 

28.305 

88.009 

75.1 

2.2 

7 

08/05/04 

19:09:26.9 

29.775 

86.985 

10.8 

2.3 

8 

08/13/04 

17:10:42.5 

29.124 

87.482 

4.0 

2.8 

9 

08/16/04 

10:20:54.8 

29.535 

85.442 

0.8 

2.3 

10 

08/16/04 

13:46:55.2 

29.528 

85.437 

10.6 

2.1 

11 

09/04/04 

17:36:54.2 

30.100 

85.580 

7.0 

3.2 

12 

09/15/04 

04:36:26.6 

30.092 

85.581 

10.5 

3.1 

13 

10/22/04 

03:28:42.5 

28.722 

87.381 

0.9 

2.7 

14 

12/08/04 

22:48:37.3 

29.796 

86.987 

0.3 

3.0 

15 

12/08/04 

22:57:20.1 

29.799 

86.984 

3.9 

2.4 

16 

12/15/04 

00:03:07.0 

28.387 

88.067 

72.1 

2.6 

17& 

01/01/05 

22:46:22.4 

29.927 

88.349 

8.2 

3.4 

18& 

01/04/05 

00:30:48.9 

29.950 

88.349 

9.7 

2.3 

19& 

01/08/05 

12:51:09.7 

29.929 

88.350 

9.4 

3.0 

20 

02/16/05 

02:33:55.6 

30.953 

85.150 

5.1 

2.6 

21& 

02/20/05 

23:41:13.0 

29.643 

88.219 

0.0 

2.2 

22& 

02/22/05 

23:38:26.5 

29.630 

88.224 

0.0 

2.5 

23& 

02/23/05 

02:03:45.6 

29.640 

88.228 

0.6 

1.9 

24& 

02/23/05 

03:30:09.2 

29.648 

88.241 

0.7 

2.6 

25& 

02/23/05 

19:53:04.1 

29.639 

88.222 

0.0 

2.2 

26 

03/19/05 

09:30:47.0 

29.945 

86.964 

7.0 

2.2 

27& 

04/30/05 

03:43:10.7 

29.639 

88.193 

0.3 

2.1 

28 

05/03/05 

08:51:23.6 

28.266 

87.737 

80.8 

3.0 

29 

05/10/05 

21:47:53.2 

28.649 

87.339 

4.7 

2.5 

30 

05/11/05 

00:45:35.5 

30.672 

85.235 

8.7 

2.8 

31 

05/11/05 

01:51:48.7 

30.667 

85.228 

12.2 

2.1 

32 

05/11/05 

04:13:00.2 

30.665 

85.232 

5.8 

4.7 

33 

05/11/05 

06:37:07.2 

28.648 

87.338 

8.5 

2.8 

34 

05/11/05 

08:25:44.6 

30.673 

85.243 

7.1 

2.5 

35& 

06/04/05 

20:44:15.6 

29.706 

88.314 

0.5 

2.2 

36 

06/09/05 

18:53:58.8 

30.125 

85.684 

0.7 

2.2 

37 

06/16/05 

18:34:32.1 

28.295 

87.648 

65.2 

2.0 

38 

07/06/05 

03:11:02.7 

28.388 

87.512 

69.6 

3.6 

39 

07/28/05 

00:21:37.9 

30.076 

85.609 

10.6 

2.4 

40 

07/29/05 

05:06:03.0 

30.674 

85.239 

12.6 

2.1 

41 

08/15/05 

10:13:46.2 

30.054 

85.662 

8.0 

2.2 
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8.  CONCLUSIONS 


•  The  Hi-CLIMB  broadband  seismic  dataset  provides  the  opportunity  for  ground  truth  location 
and  source  parameter  detennination  in  an  important  region. 

•  The  dataset  reveals  an  abundance  of  seismic  events  not  included  in  any  other  catalog.  We 
developed  modules  to  clean  up  the  automatic  event  location  database  that  eliminate  spurious 
phase  arrivals  leading  to  improved  locations,  or,  where  this  was  impossible,  remove  spurious 
events  entirely.  The  resulting  database  contains  over  22,500  events  with  over  7,900  based  on  25 
or  more  P  and  S  arrivals.  Manual  locations  of  a  small  subset  of  events  in  southern  Tibet  reveal  a 
non-systematic,  few-kilometer-scale  difference  between  locations  based  on  automatic  and 
manual  phase  picks.  Applying  additional  requirements  on  the  station  distribution  results  in  a 
high-quality  subset  of  events  with  resolved  hypocenters  that  show  a  bimodal  depth  distribution 
with  events  occurring  shallower  than  20  km  and  at  the  Moho  and  deeper. 

•  Initial  ML  detennination  indicated  systematic  underestimation  and  frequent  contamination  of 
smaller  events  by  later,  larger  events  due  to  the  high  level  of  seismicity.  Instead  of  determining 
Ml  from  S  phases,  we  implemented  a  local  magnitude  based  on  P  wave  amplitudes  ML(P)  and 
use  a  simple,  physics-based  conversion  to  obtain  regular  ML. 

•  S-picker.  During  the  last  reporting  period,  our  focus  was  on  the  implementation  of  an  automatic 
S-phase  picker,  which  uses  polarization  characteristics  specifically  geared  towards  detecting 
shear  wave  signals  and  perfonns  superior  to  a  standard  STA/LTA  detector  in  tenns  of  number  of 
detections  and  pick  accuracy.  After  optimizing  picking  parameters  and  threshold  settings  to  our 
dataset,  we  applied  the  polarization  S  picker  to  all  Hi-CLIMB  events  of  magnitude  greater  than 
2. 

•  Wavefonn  modeling  efforts  resulted  in  a  143  event  moment  tensor  database. 

•  We  identified  54  GT5-candidate  events  with  ML(P)>=2.0  of  which  50  satisfy  the  GT5  criteria 
at  the  95%  confidence  level. 

•  Two  studies  that  are  close  to  finishing  will  lead  to  further  improvements  in  event  location  and 
characterization  through  tomographic  imaging  and  attenuation  analysis. 

•  The  GT5  locations  and  moment  tensor  depths  will  contribute  to  SLBM  tomographic  efforts, 
allow  evaluation  of  depth  for  crustal  models,  and  enhance  model  accuracy  throughout  central  and 
southern  Asia  and,  in  general,  will  contribute  to  the  National  Nuclear  Security  Administration 
(NNSA)  Knowledge  Base. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


AFRL 

Air  Force  Research  Laboratory 

GSN 

Global  Seismic  Network 

GT 

Ground  truth 

NEIC 

National  Earthquake  Infonnation  Center 

PXR 

Pumqu-Xianza  Rift 

USGS 

U.S.  Geological  Survey 
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